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Abstract 

In this paper, we perform a complete analysis of the kinetic behavior of the gen- 
eral modifier mechanism of Botts and Morales in both equilibrium steady states and 
non-equilibrium steady states (NESS). Enlightened by the non-equilibrium theory 
of Markov chains, we introduce the net flux into discussion and acquire an expres- 
sion of product rate in NESS, which has clear biophysical significance. Up till now, 
it is a general belief that being an activator or an inhibitor is an intrinsic property 
of the modifier. However, we reveal that this traditional point of view is based on 
the equilibrium assumption. A modifier may no longer be an overall activator or 
inhibitor when the reaction system is not in equilibrium. Based on the regulation 
of enzyme activity by the modifier concentration, we classify the kinetic behavior 
of the modifier into three categories, which are named hyperbolic behavior, bell- 
shaped behavior, and switching behavior, respectively. We show that the switching 
phenomenon, in which a modifier may convert between an activator and an inhibitor 
when the modifier concentration varies, occurs only in NESS. Effects of drugs on the 
Pgp ATPase activity, where drugs may convert from activators to inhibitors with the 
increase of the drug concentration, are taken as a typical example to demonstrate 
the occurrence of the switching phenomenon. 
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1 Introduction 

Modifiers or effectors, ligands that bind to enzymes and thereby alter their catalytic 
activity, play a crucial role in the study of biochemical problems, e.g., enzymatic catalysis 
and metabolic pathways [l|-[5f] . Moreover, they have wide applications in pharmacology, 
toxicology, industry and agriculture. Activators and inhibitors are defined as modifiers 
that strengthen or weaken, respectively, the enzyme activity of the reaction system f^Q]. 
The enzyme activity is generally characterized in terms of the rate of product formation 
of the enzyme-catalyzed reaction in the steady state. 

Most enzyme mechanisms that involve a modifier reversibly acting on Michaelis- 
type enzymes can be regarded as a particular case of the general modifier mechanism of 
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Botts and Morales, as is depicted in Fig. Q][8]. Many theoretical biologists have studied 
the steady state and transient phase kinetics of the general modifier mechanism 



and its particular cases, in which modifiers act on Michaelis-type enzymes as competitive 



inhibitors, uncompetitive inhibitors or pure non-competitive inhibitors [l|, ll7Hl9l] 
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Figure 1: 



Segel and Martin [y, [9(] reported a steady state rate equation that is second degree 
in both substrate concentration [S] and modifier concentration [R\. They also found 
several conditions under which the rate equation can be reduced to one that is first 
degree in [S] and in [R]. Fontes et al. [7] discussed the behavior of the modifier with the 
change of substrate concentration [S] under the assumption of rapid equilibrium. Laidler 
[l9| studied the behavior of the modifier with the change of modifier concentration 
[R] under some simplifying assumptions. He also suggested definitions of competitive, 
uncompetitive and noncompetitive activation, by analogy with the generally accepted 
definitions for inhibition. 

The major differences among the contributions of these authors are the set of 
simplifying assumptions made about the steady state reached by the enzyme-catalyzed 
reaction system 20]. However, to date, there is still a lack of a complete analysis about 



the steady state kinetics of the general modifier mechanism of Botts and Morales without 
any simplifying assumptions. The major difficulty lies in the fact that the traditional 
approaches in solving this problem obscure the essence, which can hardly be realized 
without the concept of the net flux introduced in this article, of the enzyme system to 
some extent. 

It is illuminating to point out that the modifier- and substrate-binding steps are not 
dead-end reactions in the enzyme system, and so they are not necessarily in equilibrium 

We have good reasons to believe that biochemical systems in living cells generally 
operate in a state far from equilibrium. Whether the cyclic reaction mechanism in Fig{T] 
satisfies the equilibrium assumption, i.e. detailed balance, depends on whether the 
system is closed or open. A closed system will finally approach an equilibrium steady 
state, whereas an open system, driven by an external source of energy, tend to reach a 
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non-equilibrium steady state (NESS) [2ll424|. 



In this article, we remove the equilibrium assumption and provide a general analysis 
of the kinetic behavior of the modifier in NESS. Enlightened by the circulation theory of 
Markov chains 25|], we introduce the net flux into discussion and acquire an expression 



of product rate in NESS, which has clear biophysical significance. The essence of the 
general modifier kinetics is then revealed to be the competition between the equilibrium 
and non-equilibrium effects. 

So far, it is a general belief that a modifier acts as either an activator or an inhibitor 
for all its possible concentration values [R] when the substrate concentration [S] is fixed. 
However, we find that a modifier cannot be regarded as an overall activator or inhibitor 
when the reaction system is in NESS. According to our results, a particular modifier 
may convert from an activator to an inhibitor or vice versa with the change of [R]. 
More specifically, we classify the kinetic behavior of the modifier into three categories, 
which are named hyperbolic behavior, bell-shaped behavior, and switching behavior, 
respectively. The latter two kinds of behavior will never occur in an equilibrium steady 
state. 

Incidentally, drugs are typical modifiers in pharmacology. The presence of the drug 
can activate or inhibit the enzyme activity. Experimental data show that a drug can 
always act as an activator regardless of its concentration, or first act as an activator then, 
from a certain concentration value, transit to be an inhibitor 16]. Here the occurrence 



of switching phenomenon provides strong support for the argument presented in this 
paper. 



2 Methods 



2.1 Catalytic cycle 

In this article, the symbols E, S, R and P stand for the enzyme molecule, the 
substrate molecule, the modifier molecule, and the product molecule, respectively, while 
the composite symbols ES, ER and ERS represent the corresponding complexes. If 
there is only one enzyme molecule, it may convert among four states: the free (unbound) 
enzyme E, the complex ES, the complex ERS and the complex ER. Then from the 
perspective of the single enzyme molecule, the kinetics are stochastic and cyclic, as is 
shown in Fig. [2j with pseudo-first-order rate constants ai[S], a,2[R], a,s[S], and 04 [R], 
and first-order rate constants b\ = b\ + c%, 62 , &3 = 63 + C3 , and 64 . We add the rate 
constants b\ and c\, since there are two ways of transition from the complex ES to the 
free enzyme E. The rate constants 6 3 and C3 are added for the same reason. 
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Figure 2: 



2.2 Net flux 



Based on the law of mass action, we have the following kinetics equations: 
d[E] 



dt 
d[ES] 

dt 
d[ERS] 

dt 

d[E] 
dt 



-(01 [5] + a 4 [R])[E] + 61 [£75] + 64 
oi[5][£7] - (61 + o 2 [J?])[£75] + 6 2 [£7i?5], 
a 2 [i?][£75] - (62 + b 3 )[ERS] + a 3 [5] [ER] , 
a A [R][E] + 6 3 [£«5] - (64 + a 3 [5] )[£#], 



(1) 



where 61 = b[ + c\ , 63 = 63 + C3 . The above four equations constitute a system of linear 
equations with coefficient matrix 
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(2) 



Let Eq = [E] + [£75] + [ER] + [ERS] be the total enzyme concentration. Then the 
quantities \i E = [E]/E , fj, ES = [ES]/E , hers = [ERS]/E , and her = [ER]/E 
represent the probability distribution of the four enzyme states, respectively. The whole 
setup is nothing but a continuous-time Markov chain with transition density matrix Q t 
(the transpose of the matrix Q). 

It should be indicated that the enzyme-modifier and enzyme-substrate interactions 
often involve rapid binding steps followed by a slow conformational change or chemical 
step [26]. Thus, the quasi-steady approximation can be applied based on the difference 
in timescales between the catalytic cycle kinetics and the overall rate of change of bio- 
chemical reactions 27(. Assuming that the cycle kinetics represented in Fig. [2] are rapid 
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and maintain the enzyme and the complexes in a rapid quasi-steady state, we can obtain 
the steady state rate, v, of product formation for the general modifier kinetics: 



V = dt = ° l[ES] + Cs[ERS] = E °(°^es + c 3f i ERS ). (3) 
For simplicity, we normalize the product rate v and write 

V = CifJ.ES + C 3 fJ, E RS (4) 

hereinafter. 

In order to obtain the steady state concentrations [E] , [ES] , [ER] , and [ERS] , we 
set d[E]/dt = d[ES]/dt = d[ERS]/dt = d[ER]/dt = 0. Then Eq. © reduces to the 
following compact form: 

Qfi = 0, (5) 

where \i = (he, Hes, Hers, Her) 1 is a column vector representing the steady state prob- 
ability distribution of the four enzyme states. 



In the spirit of the circulation theory of Markov chains [251 ]. we introduce the 
net fluxes, namely the differences between fluxes in clockwise direction and fluxes in 
counter-clockwise direction along adjacent states of the catalytic cycle depicted in Fig. [2j 
When the reaction system reaches a steady state, the net fluxes between all adjacent 
states along the catalytic cycle are identical. Denote by J the net flux along the cycle 
E -)■ ES — >• ERS -> ER^ E. With the net flux J introduced above, Eq. © can be 
simplified to the following circulation equations: 

J = ai[S]nE ~ hHES, 

J = a 2 [R]HES ~ foHERS, ., 

(6) 

J = hHERS - a3[S\HER, 
_ J = b^HER - a 4 [R]HE- 



2.3 Equilibrium and non-equilibrium steady states 

From the perspective of a single enzyme molecule, the enzyme system depicted in 
Fig. [His in detailed balance if the net flux J along the catalytic cycle depicted in Fig. [2] 
is zero, i.e. 

ai[S]nE = hfiES, 
a.2[R]HES = b 2 HERSi 
hHERS = a 3 [S]HER, 
biHER = cl a [R]he- 

In this article, a steady state that satisfy the condition of detailed balance described 
in Eq.([7D is referred to as an equilibrium steady state. Unless otherwise stated, the 
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word 'equilibrium' appearing in this article will be always understood as an equilibrium 
steady state. However, if detailed balance is violated, then the enzyme system depicted 
in Fig. [T] is an open system that approaches a non-equilibrium steady state (NESS) 
23,0]. This is the scenario in enzyme kinetics 0, 28]. 



2.4 Steady state rate formula 

In the two extreme cases when the modifier concentration approaches zero or in- 
finity, the general modifier mechanism reduces to the single-substrate single-product 
Michaelis-Menten mechanism depicted in Fig. [3J The steady state product rate can be 
easily calculated in such cases. In fact, when the modifier concentration [R] is zero, the 
steady state rate is 

c i a i[51 / Q s 
ai[S\ +h 

Similarly, when the modifier concentration [R] is very large, the steady state rate is 

C3g3[g] 

Voo = C^ERS = rg1 , , ■ (9) 

The steady state rates Vq and Voo are collectively referred to as limiting product rates 
hereinafter. 
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Let v = her + Hers be the sum of the steady state probabilities of enzyme states 
ER and ERS. It follows from the definition of product rate (HJ), the circulation equation 
([6]), and the expressions of two limiting product rates (JSj) and © that 

, v . ,Voo VQ J 
V - Vq = (Voo - V )V + ( jy-r. (10) 

a 3 ai [S\ 

The quantity 

Av = V — Vq (11) 

in the left-hand side of Eq. ()10p is of special significance, since the sign of Av reflects the 
role of the modifier. A positive sign of Av indicates that the modifier is an activator, 
while a negative sign of Av indicates that the modifier is an inhibitor. 
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It is important to notice that the steady state rate formula ()10p has clear biophys- 
ical significance. It is easy to see from this expression that the role of the modifier is 
determined by the competition of the first term of the right-hand side which does not 
include the net flux J and the second term which includes the net flux J. 



3 Kinetic analysis 

3.1 Kinetic analysis in an equilibrium steady state 

In an equilibrium steady state, the catalytic cycle in Fig. [2] is in detailed balance 
and there is no net flux. In this case, the steady state rate formula (fTUj) suggests that 

v = v + (Voo - vq)v. (12) 

Based on the condition of detailed balance (J7J), we have the following equation in which 
the dependence of the steady state rate v on the modifier concentration [R] becomes 
extremely clear: 

VooIRI+VqK 

V= [R]+K ' (13) 

where K = (616203 + 016203 [S] )/{a\ 0263 + 010203(5]) is a constant if the substrate 
concentration [S] is fixed. 

Notice that if [S] is fixed, the product rate in equilibrium exhibits a hyperbolic 
dependence on [R\. In this case, the modifier acts as an overall activator or inhibitor 
for all its possible concentration values, depending on whether vq is smaller or greater 
than Voo. In other words, when the reaction system is in an equilibrium steady state, 
to be an activator or an inhibitor is an intrinsic property of the modifier. Fig. 2] is an 
illustration of the above conclusion. 



3.2 Kinetic analysis in NESS 

In a non-equilibrium steady state, the net flux J will no longer vanish. As is 
mentioned above, the steady state rate formula (|10p has clear biophysical significance 
and the role of the modifier is determined by the competition of the equilibrium term 

Veg = (Voo ~ Vq)v (14) 

which reflects the contribution of the equilibrium kinetics and the non-equilibrium term 

Vneq = ( ) 7^7 (15) 

a 3 ai [S\ 

which reflects the contribution of the non-equilibrium kinetics due to the non-zero net 
flux J. 
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In the following discussion, for simplicity, we assume Voo < vq. The other case of 
foo > can be discussed in the same way. Under this assumption, the equilibrium term 
v eq is negative for all possible values of [R\. It is readily verified that the equilibrium 
term v eq , which behaves like the product rate in an equilibrium steady state, exhibits an 
approximately hyperbolic dependence on [R], as is illustrated by the red curve in Fig. [5j 
On the other hand, the expression of J in the supplementary material indicates that the 
non-equilibrium term v neq exhibits a bell shape vanishing at both zero and infinity as a 
function of [R] , as is illustrated by the green curve in Fig. [5j 



Switching Behavior 
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Under certain conditions described in the next section (Sec. l4.4j ). the non-equilibrium 
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term v neq may be positive and play a dominating role compared to the negative con- 
tribution of the equilibrium term illustrated in Fig. [5j In this case, v — vq is 
positive when [R] is small and negative when [R] is sufficiently large. The quantity Av 
will change the sign somewhere in the range of the modifier concentration [R] . In other 
words, the role played by the modifier will convert from an activator to an inhibitor, as 
is illustrated by the blue curve in Fig. 

Unexpectedly, the traditional point of view that being an activator or an inhibitor 
is an intrinsic property of the modifier does not hold in NESS. The role of the modifier 
depends greatly on the modifier concentration [R]. The new phenomenon in NESS 
described above, in which the role of the modifier converts between an activator and an 
inhibitor, is referred to as switching phenomenon throughout this paper. We will discuss 
the switching behavior of the modifier in much more detail in Sec. 14.41 

Before leaving this section, we would like to point out that the switching behavior of 
the modifier may provide an explanation for the phenomenon of the side-effects caused by 
excessive intake of drugs. In fact, drugs are typical modifiers in pharmacology. If a drug 
operates in a mechanism discussed in this article that is far from equilibrium, the above 
analysis predicts that the role of the drug may convert when the drug concentration 
varies. Specifically, when a switching modifier is a drug designed to inhibit a particular 
enzyme in vitro, the concentration of the drug must be controlled within a certain range 
to maintain the inhibition effect. An exceptionally high concentration may cause the 
dangerous side-effects in which the inhibitor can produce an activation effect. 



4 Classification of the Modifier 

4.1 Preparations 

We are now prepared to classify all possible steady state behavior of modifiers 
involved in the general modifier mechanism of Botts and Morales depicted in Fig. [TJ We 
will give a quantitative criterion for various kinds of non-equilibrium kinetics in terms 
of the rate constants, Oj and 6j, and the limiting product rates, vq and Voq. To make our 
discussion friendly to those unfamiliar with tedious mathematical tools, we would like 
to present the results here and put the mathematical treatment in the supplementary 
material. 

The following three quantities play a key role in our classification: 

Al = Voo - V , 

A 2 = q (voo - v ) + k( ), f 16 \ 

o 3 ai ^ xu ; 

A 3 = (q~ q')(voo ~ v ) - k(— - — ), 

«3 «i 

where k = 01026364 — 61620304 and q,q' (q' < q) are positive parameters depending on 
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the rate constants and the substrate concentration [S]. 

Virtually, in terms of the sign of the above three discriminants, we can classify the 
behavior of the modifier with completeness and clarity. 

4.2 Hyperbolic behavior: A 2 A 3 > 

If A2A3 > 0, the steady state product rate exhibits an approximately hyperbolic 
dependence on [R\. In this case, the modifier behaves in the same way as it does in an 
equilibrium steady state except that the modification is stronger or weaker due to the 
contribution of the net flux J. The blue curves in Fig. EtA) and Fig. E^B) illustrate 
the hyperbolic behavior of the modifier as an overall activator and an overall inhibitor, 
respectively. 

4.3 Bell-shaped behavior: A1A3 < 

If A1A3 < 0, the steady state product rate exhibits a bell-shaped dependence on 
[R]. The bell-shaped behavior differs from the hyperbolic behavior in that the product 
rate will reach an extreme value which excels the limiting rate Voq with the increase of 
[R\. The red curves in Fig. E][A) and Fig.[6jB) illustrate the bell-shaped behavior of the 
modifier as an overall activator and an overall inhibitor, respectively. 

Although both the hyperbolic-behaved and bell-shaped-behaved modifiers are over- 
all activators or inhibitors for all possible values of [R] when [S] is fixed, there exist cru- 
cial differences between them. For the hyperbolic kinetics, an inordinately large increase 
of [R] is necessary to bring about even a comparatively modest change in the enzyme 
activity. For the bell-shaped kinetics, on the contrary, a modest increase of [R] will make 
the enzyme activity exceed its limit value Voq. Therefore, if the requirement for effec- 
tive regulation (activation or inhibition) of the living body is needed, then bell-shaped 
modification will be a good choice. 

4.4 Switching behavior: AiA 2 < 

If A1A2 < 0, the quantity Av will change the sign somewhere in the range of the 
modifier concentration. In this case, the role played by the modifier will convert from an 
activator to an inhibitor or vice versa. Therefore in NESS, as is pointed out in Sec. 13.21 
the effect of activation or inhibition should not be viewed as an intrinsic property of the 
modifier. 

The switching behavior of the modifier is illustrated by the green curves in Fig.[6^A) 
and Fig.[6^B). Fig. [6]( A) illustrates the transition from an inhibitor to an activator, while 
Fig. El^B) illustrates the transition from an activator to an inhibitor. 
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4.5 Classification of the modifier in terms of [S] 

Incidentally, the behavior of the modifier can be classified according to the regula- 
tion of enzyme activity by the substrate concentration [S]. In fact, the steady state rate 
v may exhibit an approximately hyperbolic or a bell-shaped dependence on [S]. The 
bell-shaped dependence differs from the hyperbolic one in that the steady state rate will 
reach an extreme value which excels the saturated rate with the increase of [S]. This 
conclusion is illustrated in Fig. \7\ where the blue and red curves represent the above 
two cases, respectively. 
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Figure 7: 
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5 Discussion 



5.1 Conditions for hyperbolic behavior 

In the following discussion, the bell-shaped behavior and the switching behavior 
of the modifier are collectively referred to as non-hyperbolic behavior. As is analyzed 
above, non-hyperbolic behavior never appears in an equilibrium steady state. In this 
section, we will discuss another case in which non-hyperbolic behavior never occurs. 

To this end, we rewrite the steady state rate formula (fT0|) as 



Av = ( Voo - vo)(u + -^) + - -)J- (17) 

a 3 [S\ [S\ a 3 ai 

It is readily verified that the first term on the right-hand side of this equation exhibits a 
hyperbolic dependence on [R\. Hence non-hyperbolic behavior will never occur as long 
as the second term vanishes. 

In an equilibrium steady state, the second term vanishes due to the zero net flux. 
Another crucial situation where the second term will vanish arises when a\ = «3, which 
implies that whether the modifier is bound to the enzyme or not will not affect the 
binding affinity of the substrate binding site. The above argument leads to the fact 
that when the binding rate of the substrate is irrelevant of the modifier, only hyperbolic 
behavior will occur, even if the enzyme system is in NESS. 

5.2 Examples of switching behavior 

As is pointed out in the last paragraph of Sec. 13.21 drugs are typical modifiers in 
pharmacology. If a drug operates in a mechanism discussed in this article that is far 
from equilibrium, the above analysis predicts that the role of the drug may convert when 
the drug concentration varies. 

ATPase activity associated with P-glycoprotein (Pgp) is characterized by three 
drug-dependent phases: basal (no drug), drug-activated, and drug-inhibited. The com- 
munication between the drug-binding site and the ATP hydrolytic site on a Pgp enzyme 
makes the reaction system a general modifier mechanism of Botts and Morales where 



ATP acts as the substrate and the drug acts as the modifier [id. l29l|. 

Experimental data show that both hyperbolic activation and switching phenomena 
occur under the experimental condition of pH 7.4 and 37° C. The four curves in Fig. [81 



which was first generated by Al-Shawi et al. [id. l3Clt] . represent the variation trends of 
enzyme activity versus concentrations of different types of drugs, namely valinomycin, 
verapamil, SL-verapamil, and colchicine, respectively. The occurrence of switching phe- 
nomenon in this experiment strongly supports the analysis of non-equilibrium kinetics 
of the general modifier mechanism presented in this paper. 



12 



500 I r 




[Drug] (M) 

Figure 8: 



5.3 Applications 



We have seen that the core of the general modifier mechanism of Botts and Morales 
is the four-state cyclic reaction system depicted in Fig. [2j Actually, the four-state 
cyclic topology is fundamental, since it models almost all possible reaction mechanism 
of proteins with two binding sides. Particularly, in living cells, biochemical processes 
that involve a four-state cyclic reaction system are very common. For example, the 
fundamental phosphorylation-dephosphorylation cycle illustrated in Fig. [9] constitutes a 
four-state cyclic reaction system [27|. 



s 



<2l[ATP] 



o,[PI], 



a 2 



KS 
'A 



6 2 [ADP] 



S*P 



s 



* 



Figure 9: 



Furthermore, in many biochemical reaction systems, proteins possess three or more 
binding sides. An eight-state cubic reaction mechanism is proposed to describe such 
kind of systems. A typical example of the eight-state cubic topology is the subunit of 
the inositol 1,4,5-trisphosphate receptor (IP3R) which is a channel located in the endo- 
plasmic reticulum that releases Ca 2+ ions. Structurally, the IP3R is a large homomeric 
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tetramer of four subunits forming a single ion-conducting channel 3l|. The gating of 
IP3R channels requires that three or all of the four subunits are at the open state 32]. 
Each subunit has one binding site for IP3 and two binding sites for Ca 2+ . Thus there 
are eight possible states for the subunit, which is illustrated in Fig. [TUJ Binding with 
IP3 'potentiates' the subunit. The two calcium binding sites activate and inactivate the 
subunit, and a subunit is activated when IP3 and the activating calcium site are bound 
but the inactivating site is unbound [33]. 



6 : open state 7 




Inhibitory Calcium 



Figure 10: 



Besides, a great number of proenzymes activation mechanisms have one or two, 
identical or different modifiers acting on the activating and the activated enzymes. It 
turns out that such activation mechanisms may be considered, really or formally, as 
particular cases of the general eight-state cubic reaction mechanism 34]. 

Mathematically, it can be proved by the theory of Markov chains that if the binding 
of one site is irrelevant to that of the other two sites, the eight-state cubic reaction 
system depicted in Fig. [10] can be reduced to the 'combination' of a uni-unimolecular 
reaction system and a four-state cyclic reaction system. Thus the four-state catalytic 
cycle discussed in this article plays a basic role in analyzing the eight-state cubic reaction 
systems. 



6 Conclusion 

In our present work, we introduce the net flux J and propose a method in dealing 
with the cyclic reaction systems. Early approaches on the general modifier mechanism 
always give so many unreadable expressions and complicated situations that the 
hidden biophysical meanings are extremely hard to clarify. In this paper, however, 
with the notion of the net flux, the expression of steady state product rate has clear 
biophysical significance. We would like to emphasize that if the reaction kinetics is 
not too complex, the approach adopted in this paper gains an advantage of being neat 
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and concise over the traditional approaches, e.g., the King-Atman method and 
numerical computation. 

In addition, we would like to point out that the single-molecule enzyme kinetics 
is now attracting broad scientific attention 35J]. The general discussion in this article 
is based on the four-state catalytic cycle depicted in Fig. [2j What we concentrate on 
is the transition among various states of a single enzyme molecule. The mathematical 
theory of non-equilibrium Markov chains and the concept of net flux along a cycle 
render possible a deep understanding of a number of interesting phenomena that never 
appear in equilibrium systems, such as the switching behavior of the modifier involved 
in the general modifier mechanism of Botts and Morales. Surely, applications of the 
non-equilibrium theory of Markov chains are not limited to the model discussed in this 
paper. It is expected to become a fundamental and useful tool in studying the single- 
molecule enzyme systems involved in signal transduction, gene regulation and synthetic 
biology. 
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Figure captions 

Fig. [TJ General modifier mechanism of Botts and Morales. The symbols E, S, R 
and P stand for the enzyme molecule, the substrate molecule, the modifier molecule 
and the product molecule, respectively, while the composite symbols ES, ER and ERS 
represent the corresponding complexes. The symbols aj, b[ and q are rate constants. 

Fig. [2j Catalytic cycle of the general modifier mechanism depicted in Fig. [TJ Since 
both reaction rates and c\ (6 3 and C3) in Fig. [JJ contribute to the transition from ES 
(ERS) to E (ER), we use the symbol 61 (63) to represent b[ + c\ (b' 3 + C3). 

Fig. El Two extreme cases of the general modifier mechanism depicted in Fig. [TJ 
When the modifier concentration approaches zero or infinity, the general modifier mech- 
anism reduces to the single-substrate single-product Michaelis-Menten mechanism de- 
picted in this figure. 

Fig.SJ Product rate (enzyme activity) in equilibrium versus modifier concentration. 
The initial enzyme activity is scaled to 1=100%. In an equilibrium steady state, the 
modifier acts as an overall activator (voo > vq) or inhibitor (voo < vq) for all possible 
values of [R], as is shown by the blue and red curves, respectively. The rate constants 
are chosen as a\ = l,a 2 = 2,63 = 2,64 = 1, 61 = 2,b 2 = 1, «3 = 2,a^ = 1, c\ = 
1.5, C3 = 1.5, [S] = 1 for the case of activation. The rate constants are chosen as 
a\ = 1, 02 = 2, 63 = 2, 64 = 1, b\ = 2, 62 = 1> «3 = 2, 04 = 1, c\ = 1.5, C3 = 0.5, [S] = 1 for 
the case of inhibition. 

Fig. EJ The switching behavior of the modifier in NESS when < vo. The red and 
green curves represent the equilibrium and non-equilibrium terms respectively and the 
blue curve represents the quantity At; = v — vq. The rate constants are chosen as a\ = 
1, a 2 = 1, 63 = 1, b 4 = 1, bi = 0.2, b 2 = 0.2, a 3 = 0.1, a 4 = 0.5, a = 1, c 3 = 4.5833, [S] = 1 
for this graph. 

Fig. [6j Enzyme activity in NESS versus modifier concentration. (A) Enzyme 
activity in NESS versus modifier concentration when > vq. The initial enzyme 
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activity is scaled to 1=100%. The rate constants are chosen as a\ = 20, a 2 = 20,63 = 
10,64 = 10,6i = 1,62 = 1,03 = 50,a4 = 10, ci = 1,03 = 2.3, [S] = 1 for the hyperbolic 
behavior. The rate constants are chosen as a\ = 1,02 = 10,63 = 20,64 = 30, 61 = 
1,62 = 1,«3 = 10,a4 = l,ci = 1, C3 = 21, [S] = 1 for the bell-shaped behavior. The rate 
constants are chosen as a\ = 1,0,% = 10,63 = 10,64 = 10, 61 = 1,62 = 1,03 = 0.5, 04 = 
1, ci = 1,03 = 3, [S] = 1 for the switching behavior. (B) Enzyme activity in NESS 
versus modifier concentration when Voo < vq. The rate constants are chosen as a\ = 
10, a 2 = 10, 63 = 10, 64 = 10, 61 = 1, 6 2 = 1, a 3 = 10, a 4 = 10, c x = 1, c 3 = 0.9, [S] = 1 for 
the hyperbolic behavior. The rate constants are chosen as a\ = 1, a 2 = 10, 63 = 50, 64 = 
50, 61 = 1, 62 = 1, a 3 = 20, 04 = 1, ci = 1, C3 = 0.9, [S] = 1 for the bell-shaped behavior. 
The rate constants are chosen as a\ = 1,02 = 1,63 = 1,64 = l,6i = 0.2, 62 = 0.2, 03 = 
0.1, 04 = 0.5, ci = 1, C3 = 4.6, [S] = 1 for the switching behavior. 

Fig. [71 Product rate (enzyme activity) in NESS versus substrate concentration. 
The saturated enzyme activity is scaled to 1=100%. The blue and red curves represent 
the hyperbolic and bell-shaped dependence, respectively. The rate constants are chosen 
as ai = 1, 02 = 1 ; 63 = 1, 64 = 1, 61 = 1, 62 = 1, 03 = 1, 04 = 2, ci = 390, C3 = 10, [R] = 1 
for the hyperbolic dependence. The rate constants are chosen as a\ = l,a 2 = 1,63 = 
1,6 4 = 1,61 = 0.2,62 = 0.2, a 3 = 0.1, a 4 = 0.2, d = 1150, c 3 = 10, [R] = 1 for the 
bell-shaped dependence. 

Fig. [3 Effects of drugs on Pgp ATPase activity. The ATPase activities of MDR1 
reconstituted into mixed lipid vesicles were measured as a function of added drug concen- 
tration at pH 7.4 and 37° C 0,0 ■ black circles, valinomycin; inverted black triangles, 
verapamil; gray diamonds, SL-verapamil; gray triangles, colchicine. Data for verapamil, 
SL-verapamil, and colchicine are from Omote and Al-Shawi 30]. 

Fig. [9l A typical cellular biochemical switch consisting of a phosphorylation- 
dephosphorylation cycle. The substrate molecule S may be a protein or other signaling 
molecule. If S is a protein then the phosphorylation of S is catalyzed by a protein kinase 
K and the dephosphorylation is catalyzed by a protein phosphatase P. The entire cycle 
is accompanied by the reaction ATP =^ ADP+PI. 

Fig. d0j Transition diagram for the eight-state subunit of IP3 receptor channel. 
When the system is at state 6 (one IP3 and one activating Ca 2+ bound), the subunit is 
activated. 
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